A search for radioactive ^^Al in the nova-Uke variable V4332 
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- - ' ABSTRACT 

We have searched for the important radioactive isotope ^®A1 in the nova-Hke source V4332 
Sgr. Recent results from gamma ray astronomy show that there is pervasive emission of the 
1.809 MeV gamma ray photon, arising from the decay of ^^Al to ^^Mg, from aU over the galactic 
plane. Though the sites from where this emission originates are not clearly established, novae 
are believed to be an important contributing source. In this context, V4332 Sgr presented a rare 
opportunity to observationally investigate whether novae or novae-like sources synthesize ^^Al 
and to what extent. Strong AlO bands in the near-IR have been reported in this object recently. 
As molecular bands of different isotopic compositions are readily resolved spectroscopically (e.g. 
^^CO and ^^CO), it was thought that the components of AlO associated with ^^Al and stable ^^Al 
could be detected as separate bands. Our spectra indicate that there is no strong presence of ^^Al 
in V4332 Sgr. A reliable upper limit of 0.10 for the ^^Al/^^Al ratio is determined which constitutes 
the first observational constraint for this ratio in a potential ^^Al producing source. While V4332 
Sgr is not a typical nova, its outburst amplitude and light-curve behaviour bear close similarity 
to that of novae. Hence, although the results from V4332 Sgr cannot be directly extended to 
novae in general, the limit on the observed ^^Al/^^Al ratio could be a useful input in constraining 
rather uncertain nucleosynthesis models for the production of ^^Al in novae/novae-like sources. 
By comparing the observed ^^Al/^^Al ratio in y4332 Sgr with that expected in classical novae 
it appears unlikely that the progenitor of V4332 Sgr is an Oxygen-Neon-Magnesium (ONeMg) 
white dwarf. 

Subject headings: novae, cataclysmic variables - nuclear reactions, nucleosynthesis, abundances - stars: 
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1. Introduction 

^^Al is an important radionuclide for tracing 
and studying galactic radioactivity. ^® Al - the iso- 
tope of ^^Al - is unstable with a mean life of T26 
= 1.05x10^ years and decays to ^^Mg emitting a 
gamma ray photon of 1.809 MeV. It is through the 
detection of this 1.8 MeV 7 ray photon that ^^Al 
manifests itself. The discovery of ^^Al towards 
the galactic center - deemed to be epochal in na- 
ture - was made by the High Energy Astronomical 
Observatory (HEAO-C) in 1979-80 (Mahoney et 
al. 1982) and subsequently confirmed by a few 
other missions (Prantzos & Diehl 1996). How- 
ever, these missions were limited by inadequate 
spatial resolution ( typically a few tens of degrees) 
and could not therefore locate the precise sites of 
the 1.8 MeV 7 ray emission. A major advance 
in mapping the 1.8 MeV 7 ray emission from the 
galaxy emerged from recent studies by the Comp- 
ton Gamma ray Observatory (CGRO) (Diehl et 
al. 1995). The CGRO mapping at 1.8 MeV - ob- 
tained at a vastly improved resolution (telescope 
angular response is ^ 3.8 degrees) but which is still 
below that required for a precise point source asso- 
ciation - traces ^^Al emission in localized regions 
lying mostly along the galactic plane. A diffuse 
and irregular emission along the galactic plane is 
observed favoring a massive star origin. Super- 
novae, novae. Wolf Rayet stars and AGB (Asymp- 
totic Giant Branch) stars are considered the prime 
candidates for the origin of the 1.8 MeV emission 
but this has to be established more firmly. Hence 
a convincing detection of ^^ Al in any of these class 
of objects would give more definite evidence on the 
exact sites of ^^Al production. Further, a detec- 
tion of ^^Al would have important implications for 
current models of nucleosynthesis. Such models 
give estimates - which are believed to be uncer- 
tain - of the ^^Al/^^Al ratio that can be expected 
in the different classes of ^^Al producing sources. 
An observational constraint on the ^®A1/^^A1 ra- 
tio would therefore be vitally important for nu- 
cleosynthesis modeling. The efforts from CGRO 
are being extended by the gamma ray satellite IN- 
TEGRAL which has high-resolution gamma ray 
spectroscopy as one of its main goals. Thus the 
study of the origin and distribution of ^^Al in 
the galaxy, though in a nascent stage, is assum- 
ing significant importance. It may also be men- 
tioned that ^®A1 plays an important role in so- 



lar system studies where its decay to ^^Mg is in- 
voked to explain the excess ^^Mg found in mete- 
oritic samples (the ^^Mg anomaly). It has been 
proposed that the ^^Al and several other short- 
lived nuclides present in the early solar system 
were freshly synthesized nuclides injected into the 
proto-solar cloud from a nearby source such as a 
nova/supernova/ AGB star ( Cameron et al. 1995; 
Arnould et al. 1997; Goswami & Vanhala 1998). 
Alternatively, the possibility of energetic-particle 
interactions - as a source of ^^Al production - has 
also been proposed (Gounelle et al. 2001). 

The possibility of a first, spatially-resolved and 
alternate method of detecting ^^Al (apart from 
detecting the 1.8 MeV line) in an astronomical 
object suggested itself when several strong A-X 
bands of the AlO radical were detected in V4332 
Sgr (Banerjee et al. 2003). By virtue of hav- 
ing undergone a nova-like outburst, V4332 Sgr 
satisfied the requisite criterion of belonging to a 
class of objects predicted to be a likely site of 
^^Al creation. It was also recognized that molec- 
ular bands, arising from different isotopic compo- 
sitions, can be easily resolved spectroscopically at 
an intermediate resolution. This is not the case 
for atomic lines where the isotope shifts are too 
small to detect. A frequently encountered exam- 
ple, of isotopically differing molecular bands be- 
ing resolved, arises in the first overtone bands of 
carbon monoxide at 2.3/im. Here, the signatures 
of the ^^CO and ^'^CO bands are routinely seen 
to be distinctly separate in the spectra of cool 
stars. We thus felt that the two components of 
AlO viz. the radioactive component ^^AlO (in 
case it is present at detectable limits) and the sta- 
ble component ^^AlO could be expected to be dis- 
tinctly resolved in V4332 Sgr. Considering that 
Al bearing molecules are seen very rarely in as- 
tronomical sources, V4332 Sgr presented a unique 
chance to detect the presence of ^®A1. Further, 
since the AlO bands in V4332 Sgr are strong, even 
a low strength of ^^Al vis-a-vis ^^Al in the object 
would still permit the ^^AlO bands to be detected. 
The detected strength of the putative ^^AlO bands 
would help to establish the ^®A1/^^A1 ratio in the 
source. While all the A-X bands of AlO are ex- 
pected to exhibit the ^^AlO component, the (2,0) 
band was chosen for the present study because it is 
strong and the adjacent spectral region around the 
wings of the band - where the ^^AlO signature is 



expected - is relatively uncontaminated with other 
features. Further, the wavelength separation be- 
tween the ^^AlO and ^^AlO components of the 
(2,0) band is larger than that for many of the other 
bands and therefore expected to be more easily re- 
solved. Specifically, the expected band-origin for 
the A'^113/2 components for ^"^AlO and ^^AlO are 
at 1.50351/im and 1.50155/im respectively ( a sep- 
aration of 19. 6A); for the A^Ili/2 component the 
positions are at 1.47492/im and 1.47325/im ( a sep- 
aration of 16. 7A) respectively. Such separations 
can be resolved at a spectral resolution of 1000 or 
higher (we have used a resolution of 3800 in our 
observations) . 

It was also felt that the present study could 
give additional insights into the outburst mecha- 
nism for V4322 Sgr. V4332 Sgr is intrinsically a 
very interesting source having a rich optical and 
near-IR spectrum (Banerjee et al. 2003; Baner- 
jee & Ashok 2004). It erupted in 1994 with a 
~ 9.5 magnitude brightening in the V band fol- 
lowed by a slow decline. The distance to the ob- 
ject has been estimated to be 300 pc (Martini et 
al. 1999). The puzzling aspect of the object was 
its rapid cooling to a M giant/supergiant which 
was uncharacteristic of a classical nova (Martini 
et al. 1999). Classical novae generally evolve to- 
wards the hot nebular/coronal state. A similar be- 
havior of evolving towards cool temperatures like 
V4332 Sgr, has been seen recently in the spectacu- 
lar outburst the eruptive variable V838 Mon (Mu- 
nari et al. 2002; Banerjee & Ashok, 2002a; Bond 
et al. 2003 - for the striking light echo around the 
object). It has been suggested that V838 Mon, 
V4332 Sgr and M31 RV (another red variable that 
erupted in M31; Rich et al. 1989) could belong 
to a new class of objects with outburst properties 
significantly different from those of conventional 
eruptive variables like classical, symbiotic, recur- 
rent novae or born again AGB stars (Munari et al. 
2002). There is an imperative need to understand 
the cause of the explosion in such eruptive vari- 
ables. Plausible explanations for the outburst in- 
voke merger of stars (Soker & Tylenda 2003 ) and 
a star capturing its encircling planets (Retter & 
Marom 2003). The conventional scenario for out- 
bursts in classical novae - that of a thermonuclear 
runaway on a white dwarf surface accreting matter 
from its secondary companion - does not appear 
to hold (Bond et al. 2003) but a variant of it may 



be applicable to V4332 Sgr given the uniqueness 
of the object. Since in a classical novae outburst, 
theoretical estimates indicate that the ^®A1/^^A1 
ratio can be significantly enhanced to values even 
exceeding unity , (Gehrz et al. 1999 and refer- 
ences therein; Jose & Hernanz 1998; Prantzos & 
Diehl 1996) the detected strength of ^^Al relative 
to ^^Al in V4332 Sgr could shed further light on 
the progenitor of V4332 Sgr which underwent the 
outburst. 

2. Observations 

H band spectra were obtained on 22 September 
2003 with the 3.8m United Kingdom IR Telescope 
(UKIRT), Hawaii and the 1-5 micron Imager Spec- 
trometer (UIST) using a 1024x1024 InSb array. 
Observations were done by nodding the telescope 
on two positions separated by ~ 12" along a 2 pixel 
(0.24") wide slit resulting in a spectral resolution 
of 3800. Flat fielding was done with a black body 
mounted inside UIST and spectral calibration was 
done using an Argon lamp. The standard star BS 
6998 was observed at similar air-mass as V4332 
Sgr and its spectrum was used for ratioing the ob- 
ject spectra for removal of telluric lines. Table 1 
gives the details of the spectroscopic observations. 

3. Results 

The observed spectra in the two dithered po- 
sitions are shown in Figure 1. We show these 
spectra separately to bring out the genuineness 
and repeatability of several finer features that are 
evident in the bands. The spectra here may be 
compared with the lower-resolution spectrum (R 
= 450) of the same bands as given in Banerjee at 
al. (2003) - more detail is discernible in the present 
spectra. The AlO bands described here arise from 
transitions between the A-X electronic states of 
the AlO radical (Launila & Jonsson 1994). The 
A state is a doublet with the spin components 
A^Il2/2 and A^Iii/2- The observed spectrum here 
is typical of electronic-vibration transitions in di- 
atomic molecules (Herzberg 1950). The fine struc- 
ture seen in the spectra arise from the rotational 
lines of the P, Q and R branches. The R branch 
is manifested most clearly since it gives rise to the 
prominent band-head (blueward of the band ori- 
gin and marked in Figure 1) where the R branch 
lines crowd together. The actual structure of the 



the R band-head is complex as it consists of two 
or more branches (Launila & Jonsson 1994). It is 
around this position of the R branch band-head 
that the contribution of ^^AlO bands is expected 
to be seen. Since hints of an additional component 
were seen at around this position in our earlier 
data, we felt that it could be a tentative signa- 
ture for ^^AlO. Specifically we refer to the J band 
spectrum for the yl^ni/2 component of the (4,0) 
band which shows a sloping blue wing at 1.2256/im 
(Banerjee at al. 2003; the blue wings of all the 
other AlO bands rise sharply - typical of molecu- 
lar band structure). A similar J band spectrum 
for the (4,0) band taken at a higher resolution of 
2050, but not included in Banerjee et al. (2003), 
decomposes the sloping blue wing into two distinct 
components. Thus, it was felt that the potential 
existed to try and make a more detailed search for 
26 AlO in V4332 Sgr. 

Because of the good reproducibility in the two 
spectra of Figure 1, we have averaged them to 
increase the S/N and the resultant spectrum is 
shown in Figure 2 (top panel) by the grey line. 
A computed model spectrum for AlO - assuming 
a pure ^^AlO component and no ^^AlO - for the 
(2,0) band for a rotational temperature of 200K is 
overlaid on this. As may be seen, there is a good 
match between the model and the observed data. 
The broader features, like the excellent wavelength 
matching of the R branch band- head and the over- 
all band structure are well reproduced. In ad- 
dition, many of the finer features due to the in- 
dividual rotational transitions arc well matched. 
The observed AlO spectrum of V4332 Sgr appears 
to be quite pure. However, there is a broad ob- 
served feature at 1.526/im - marked with an as- 
terisk - which does not match the model spec- 
trum and which appears to be due to some other 
species apart from AlO. The synthetic AlO spec- 
tra were created starting from the known molecu- 
lar constants of both states involved in the tran- 
sition of the experimentally known ^''AlO iso- 
topomer. Corrections were subsequently applied 
on all these constants, according to standard pro- 
cedures for isotopic corrections, and correspond- 
ing 26 AlO spectra were generated using a syn- 
thesis routine for molecular spectra (available at 
http://kurslab.physics.kth.se/ oUi/superb/). 

In the lower panel of Figure 2, we show four 
plots which indicate how the model spectrum for 



pure ^'AlO would be affected by including an 
26 AlO component. The fractional contribution 
of 26 Al in each plot is indicated in terms of the 
26A1/2^A1 ratio used to construct each curve. To 
avoid crowding, we do not show the pure 26 AlO 
spectrum - it has a structure essentially similar 
to the 2'^ AlO spectrum but laterally shifted to a 
slightly shorter wavelength (the relative displace- 
ments of the band origin are 19.6 and 16. 7A for the 
(2,0) doublets respectively). From the plots of the 
lower panel of Figure 2, it is seen that the most dis- 
tinctive signature for the presence of 26 AlO would 
be its R branch band- head labeled as feature C. 
At higher values of the 26A1/2'^A1 ratio (0.3,0.5 and 
1) it is predicted to be strong. Even at a value as 
low as 0.1, this feature is discernible in the com- 
puted plots for the A^Ii'i/2 and A^Iii/2 compo- 
nents. Hence there should have been some sign 
of it in the observed spectra if 26 Al was present 
at significant strength in V4332 Sgr . The inten- 
sity of the feature B with respect to A, also in- 
creases with increasing values of the 26A1/2'^A1 ra- 
tio. The ratio of these two features in the observed 
spectra would appear to be best replicated for a 
26A1/2'^A1 ratio of 0.1 or lower. Thus, it would ap- 
pear that 26 Al is not seen at appreciable strength 
in V4332 Sgr. To get a more rigorous estimate on 
the upper limit on the strength of 26 Al, we applied 
a Kolmogorov-Smirnov (KS) test to the observed 
and model data (for different 26Al/2'i'Al values). 
The KS test was applied to a restricted portion of 
the data between 1.47/xm to 1.48/im. This section 
of data was chosen because the strength of fea- 
tures A, B and C are expected to be considerably 
influenced by the presence of 26 Al. The KS test, 
in this region, should therefore be sensitive to the 
signature of 26 Al. From the analysis we find a high 
probability of association (P) between the model 
and observed data when 26 Al is not present (i.e. 
P = 0.91 for 26A1/27A1 = 0). As the 26A1/27A1 
ratio is increased to values of 0.05, 0.1, 0.15 and 
0.20 respectively, P falls to values of 0.80, 0.20, 
0.033 and 0.001 respectively. The computed prob- 
ability values indicate that 26 Al is definitely not 
present at large strengths of 15-20 percent. An 
upper limit of ~ 10 percent is suggested based on 
the analysis. The present data do not permit set- 
ting a lower limit on the strength of 26 Al in V4332 
Sgr. Spectra obtained at a higher resolution than 
the present studies - and at similar or greater S/N 



- where the individual rotational lines are better 
resolved could give a more definite answer. 

4. Discussion 

A positive detection of ^®A1 would have con- 
stituted definite observational confirmation of the 
theoretical prediction that novae explosions syn- 
thesize ^^Al. Even otherwise, the fairly reliable 
limit of the ^® Al/^^ Al ratio that we obtain has im- 
portant implications. The nucleosynthesis of ^^Al 
in novae and other ^^Al producing sources is a 
very complex field. It is recognized that consid- 
erable uncertainty exists in crucial areas related 
to nuclear reaction networks and nuclear reaction 
rates thereby leading to uncertainties in the pre- 
dicted value of the ^®A1/^^A1 ratio for the differ- 
ent sources (Prantzos & Diehl 1996; Gehrz et al. 
1998; Busso, Gallino & Wasserburg 1999). As 
an example we note that the reaction rates for 
the key reactions for the production and destruc- 
tion of 26 Al viz. 25Mg(p,7)26Al and 26Al(p,7)27Si 
(proton capture on ^^Mg and ^^Al) have been re- 
vised by a factor of 5 in recent times (Prantzos 
& Diehl 1996). There are additional uncertain- 
ties involving the treatment of convection in no- 
vae. An observational constraint on the ^^Al/^^Al 
ratio would therefore be extremely important for 
nucleosynthesis models. The predicted model val- 
ues for this ratio for classical novae are as follows. 
For CO novae (the outburst occurs on a Carbon- 
Oxygen white dwarf (WD) with Mwd < l-2Mo), 
the ^^Al/^^Al ratio depends strongly on the WD 
mass. For WD masses between 0.6 to 1.2M0, 
the recent results of Starrfield et al. (1997) show 
that this ratio varies between 0.01 to 0.6. For the 
ONeMg novae (which comprise ^ 30-40 percent 
of all classical novae) the WD is more massive (> 
1.2Mo) and a much greater ratio for ^^Al/^^Al be- 
tween 1-10 can be expected (Starrfield et al. 1997; 
Prantzos & Diehl 1996). The detailed nucleosyn- 
thesis models by Jose & Hcrnanz (1998) indicate 
that the ^^Al/^^Al ratio can lie between 0.2-0.4 
for ONe novae and between 0.15-0.56 for CO no- 
vae. The increased yield of ^^Al in ONeMg no- 
vae is essentially due to the existing magnesium 
rich environment in the outbursting progenitors 
in these objects (^^Mg and ^^Mg are the seed nu- 
clei for 26 Al generation). As can be seen a signifi- 
cant value of the ^^Al/^'^Al ratio can be expected 
in novae. Novae are estimated to contribute 0.1 



to O.4M0 (Jose, Hernanz & Coc 1997) of the ap- 
proximately 2M0 of 26 Al present in the galaxy (a 
value based on the HEAO-C and CGRO results). 
While V4332 Sgr is not a typical novae, and the 
predicted 26A1/2'^A1 ratio for classical novae can- 
not be strictly extrapolated here, it must still be 
noted that it underwent a single large outburst like 
a classical nova. 

It would be useful to mention the expected 
26A1/27A1 ratio in some of the other potential ^6 Al 
producing candidates. In Wolf Rayet stars this ra- 
tio is 2-5x10^2 and in type II supernovae it is 5- 
7x10^'^ . In AGB stars, it is expected to be a few 
times 10^'^ if 26 Al is produced only in the Hydro- 
gen burning shell. This ratio is uncertain, but can 
be larger than unity, in case hot bottom burning 
is also considered (Prantzos & Diehl 1996). For 
V4332 Sgr specifically, the lack of 26 Al at appre- 
ciable strength, indicates that its progenitor is un- 
likely to be a ONeMg white dwarf. As mentioned 
earlier, there is keen interest in understanding the 
nature of the progenitor and the outburst mecha- 
nism for this object and its probable analog V838 
Mon. The present results give additional insights 
in this direction. It may also be noted that 22Na 
(mean life ^ 3.75 years) is another important ra- 
dioisotope whose production is expected to be sig- 
nificantly enhanced in nova outbursts (Starrfield 
et al. 1997, Gehrz et al. 1998). Emission lines 
from the resonance doublet of Nal at 5890, 5896A 
have been detected with considerable strength in 
V4332 Sgr (Banerjee & Ashok 2004). If a signifi- 
cant fraction of the observed Nal emission is from 
22Na, then because of its short mean life, a consid- 
erable reduction in the intensity of the Nal dou- 
blet can be expected (vis-vis other emission lines 
seen in the spectrum) over a timescale of a few 
years. The object should be monitored to detect 
such changes. 22Na is also of interest for 7 ray line 
study by INTEGRAL because of the 1.275 MeV 
photon emitted by its decay to 22 Ne. 
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Fig. 1. — The two observed spectra of V4332 Sgr 
covering the (2,0) band of the A-X band system of 
AlO. The spectra were obtained from UKIRT on 
22 September 2003 at a resolution of ^ 3800. 




Wavelength (iTiicrons) 



Fig. 2. — The top panel shows the sum of the two 
spectra of Figure 1 (in gray) with a slope correc- 
tion applied. Overlaid on this is a model spectrum 
of ^^AlO (black line) for an assumed rotational 
temperature of 200K. A good match is found ex- 
cept for the broad feature marked with an aster- 
isk which is due to an unidentified species and not 
AlO. The bottom panel shows the effect on the 
pure ^^AlO spectrum of the top panel by includ- 
ing a contribution from ^® Al at different strengths. 
Please refer Section 3 for further details. 



Table 1: Log 


of spectroscopic observations 




Object 


UT Exposure Integration 
Time(s) Time(s) 


Mean 
Airmass 


V4332 Sgr 
BS 6998 


5.558 240 2880 
5.067 3 48 


1.35 
1.32 



